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The calculation of quark propagators for Ginsparg- Wilson-type Dirac operators is costly and thus limited to 
a few different sources. We present a new approach for determining spatially optimized operators for lattice 
spectroscopy of excited hadrons. Jacobi smeared quark sources with different widths are combined to construct 
hadron operators with different spatial wave functions. We study the Roper state and excited p and n mesons. 



1. Optimized quark sources 

The low-lying hadron spectrum shows a few 
features which are fingerprints of QCD. In the me- 
son sector the pion occurs as an almost-Goldstone 
boson with its squared mass vanishing propor- 
tional to the quark mass in contrast to all other 
mesons. The observed ordering of the lowest pos- 
itive, 1/2+, iV(1440), and negative parity excita- 
tions of the nucleon, l/2~, iV(1535) is 'unnatu- 
ral'. A physical picture based on linear confine- 
ment, Coulomb and color-magnetic terms, always 
arranges the first radial excitation above the first 
orbital excitation, i.e. the excited states have al- 
ternating parities. 

Whereas ground state spectroscopy on the lat- 
tice is by now a well understood physical prob- 
lem with impressive agreement with experiment, 
the lattice study of excited states is not so far 
advanced. In a lattice calculation the masses of 
excited states show up in the sub-leading expo- 
nentials of Euclidean two point functions. A di- 
rect fit of a single correlator is cumbersome since 
the signal is strongly dominated by the ground 
state. Also with methods such as constrained fits 
PP or the maximum entropy method |2] one still 
needs very high statistics for reliable results 1 1 . 
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An alternative method is the variational method 
|Jj where one diagonalizes a matrix containing all 
cross correlations of a set of several operators with 
the correct quantum numbers. For a large enough 
and properly chosen set of basis operators each 
eigenmode is then dominated by a different phys- 
ical state. After normalization the largest eigen- 
value gives the correlator of the ground state, the 
second-largest eigenvalue corresponds to the first 
excited state, and so on. 

It is important to optimize the spatial proper- 
ties of the interpolating operators. An example 
for this fact is the Roper state where the varia- 
tional method, based on nucleon operators that 
differ only in their diquark content but have the 
same spatial wave function, did not lead to suc- 
cess . It can be argued that a node in the radial 
wave function is necessary to capture reliably the 
Roper state or other radially excited hadrons. Re- 
cently w e demonstrated that an elegant solu- 
tion is to combine Jacobi smeared quark sources 
with different widths to build the hadron oper- 
ators and compute the cross-correlations in the 
variational method. We find good effective mass 
plateaus for the first and partly the second ra- 
dially excited states. The propagators are then 
fitted using standard techniques. 

Already in jS] Jacobi smeared sources 
were combined with point sources and cross- 
correlations studied in similar spirit (see also [§])• 



2 



The technique of Jacobi smearing is well known 
|8I10| . The smeared source lives in the timeslice 
t = and is constructed by iterated multiplica- 
tion with a smearing operator H on a point-like 
source. The operator H is the spatial hopping 
part of the Wilson term at timeslice 0; it is trivial 
in Dirac space and acts only on the color indices. 
This construction has two free parameters: The 
number of smearing steps N and the hopping 
parameter n. These can be used to adjust the 
profile of the source. Here we work with two 
different sources, a narrow source n and a wide 
source w with parameters given by 

n: AT = 18, k = 0.210, a/2 « 0.27 fin , , > 
w: N = 41, k = 0.191, er/2 w 0.41 fm. { ' 

N and k were chosen such that the profiles ap- 
proximate Gaussian distributions with the indi- 
cated half- widths 0. We remark that the two 
sources allow the system to build up radial wave 
functions with and without a node. The param- 
eters were chosen such that simple linear combi- 
nations c n n + c w w of the narrow and wide pro- 
file approximate the first and second radial wave 
functions of the spherical harmonic oscillator: 
The coefficients c„ ~ 0.6, c w ~ 0.4 approximate 
a Gaussian with a half- width of a/2 ~ 0.33 fm, 
while c„ ~ 2.2, c w ~ —1.2 approximate the corre- 
sponding excited wave function with one node. 

The final form of the wave function is deter- 
mined through the variational method 5 . In 
this approach one computes a complete correla- 
tion matrix of operators Oi,i — 1,2, ...R that 
create from the vacuum the state which one wants 
to analyze. The eigenvalues \( k >(t) of the cor- 
relation matrix behave as \( k ' (t) oc Mk [1 + 
0(e- tAMk )}, where AM k is the distance of M k to 
nearby energy levels. The hadron sources we use 
for the correlation matrix are constructed from 
the narrow and wide quark sources. 

2. Excited nucleon signals 

For our quenched calculation we use the chi- 
rally improved Dirac operator |11| . It is an ap- 
proximation of a solution of the Ginsparg- Wilson 
equation which governs chiral symmetry on the 
lattice. This operator is well tested in quenched 



ground state spectroscopy ^2] where pion masses 
down to 250 MeV can be reached at a consider- 
ably smaller numerical cost than needed for exact 
Ginsparg- Wilson fermions. For ground states the 
chirally improved action shows very good scaling 
behavior. The gauge configurations were gener- 
ated on a 12 3 x 24 lattice with the Liischer-Weisz 
action The inverse gauge coupling is fi = 7.9, 
giving rise to a lattice spacing of a = 0.148(2) fm 
as determined from the Sommer parameter |14|. 
The statistics of our ensemble is 100 configura- 
tions. We use 10 different quark masses m rang- 
ing from am = 0.02 to am = 0.20. 

Our analysis is based on the interpolator 
£abc(uaC"f5db)u c . Each of the three quarks can 
be smeared either narrow (n) or wide (w). This 
gives 8 possible combinations (nnn, nnw, etc.). 
From a subset of 4 of these operators (after pro- 
jection to definite parity) we calculate the corre- 
lation matrix C (t) which we then use in the vari- 
ational method. The exponential decay of three 
eigenvalues is clearly identified. We identify these 
signals with the nucleon, the Roper state and the 
next positive parity resonance AT (1710). A de- 
tailed discussion of further checks on the correct 
identification of the Roper state can be found in 
[7] (see also [3] concerning the problem of nucleon- 
rj' ghost contributions). 

3. Excited meson signals 

As another test of our approach we discuss 
the 7r- and p-mesons and their radial excitations. 
For the p we use the interpolators u(x) 7i d(x) 
and u{x) 74 7, d(x), for the pion u(x) 75 d{x) and 
u[x) 747s d[x). Again we use wide and narrow 
quark sources for both interpolating fields, corre- 
sponding to 3 operators each (the combinations 
nw and wn give identical correlators and one of 
them can be omitted). When diagonalizing the 
3x3 matrix with either interpolator we see a pro- 
nounced exponential decay only for the two larger 
(in magnitude) eigenvalues, X^\t) and \( 2 \t). 
The smallest eigenvalue X^ 3 \t) does not show a 
clear effective mass plateau. This is an indica- 
tion that this eigenvalue couples to an unphysical 
quenched ghost state |15I16I3| . The final results 
for the masses as a function of the quark mass are 
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Figure 1. Masses of p(770), p(1450), tt(140) and 
7r(1300) as a function of the quark mass (the ex- 
perimental data were converted to lattice units 
with the Sommer parameter scale). 



shown in Fig. ^ 

We find that the ground state meson masses ap- 
proach their experimental values reasonably well. 
The excited state masses are considerably above 
their experimental values. There is, however, a 
plausible reason for this behavior. The sizes of 
hadrons which are not, or only weakly, affected 
by spontaneous chiral symmetry breaking can be 
estimated from the known string tension, which 
is approximately f GeV/fm. Hence the size of the 
excited mesons should be larger than the ground 
state, about 1.5 fm. Thus the size of our lattice 
(1.8 fm) is clearly not enough for a precise mea- 
surement of e.g. the p(1450) mass. The finite size 
effect cannot be neglected for the excited state 
since it apparently shifts the measured mass up 
as compared to the experimental value. 



A crucial test of our method is to check whether 
indeed the ground state is built from a nodeless 
combination of our sources and the excited states 
do show nodes. This question can be addressed 
by analyzing the eigenvectors of the correlation 
matrix. This has been done in [7] and indeed 
confirms the expectation. 
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